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Abstract. In this paper we show, in dimension n > 3, that 
Q> [ knowledge of the Cauchy data for the Schrodinger equation in the 

presence of a magnetic potential, measured on possibly very small 
subsets of the boundary, determines uniquely the magnetic field 
p I ■ and the electric potential. We follow the general strategy of [7] 

I using a richer set of solutions to the Dirichlet problem that has 

• ' been used in previous works on this problem. 

1. Introduction 

Let Q C R" be an open bounded set with C°° boundary, we are 
' interested in the magnetic Schrodinger operator 

(1.1) CA,,ix, D) = J2iDj + Mx)f + q{x) 

O . i=i 

O' =D'^ + A-D + D-A + A^ + q 

\ with real magnetic potential A = (^j)i<j<n £ C^(r2; R*^) and bounded 

' electric potential q G L^{VL). As usual, D = —iV. In this paper, we 

C ■ always assume the dimension to be > 3. Let us introduce the 



> 

X 



Assumption 1. is not an eigenvalue of the magnetic Schrodinger 

X ■ operator Ca n ■ H^Q) n H^{n) ^ L'^^Q) . 

H 

Let u be the unit exterior normal. Under the assumption Q the 
Dirichlet problem 



:i.2) 



^A,qU = 

u\an = f eHHdQ) 

has a unique solution in H^{Q), and we can introduce the Dirichlet to 
Neumann map (DN) 

associated to the magnetic Schrodinger operator £^ ^ with magnetic 
potential defined by (jl.ip . 

1 
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The inverse problem we consider in this paper is to recover infor- 
mation about the magnetic and electric potential from the DN map 
measured on subsets of the boundary. As was noted in [I^, the DN 
map is invariant under a gauge transformation of the magnetic poten- 
tial: it ensues from the identities 

(1.3) e~**£A,ge** = £yi+v*,q, e~*"^A/A,ge*"^ = A/a+v^-,?, 

that NA,q = J^A+S7^,q when \I' G C^{Q) is such that "^lan = 0. Thus 
AfA,q carries information about the magnetic field^ B = dA. Sun showed 
in [12j that from this information one can determine the magnetic field 
and the electric potential if the magnetic potential is small in an appro- 
priate class. In [H] the smallness assumption was eliminated for smooth 
magnetic and electric potentials and for and compactly supported 
magnetic potential and L°° electrical potential. The regularity assump- 
tion on the magnetic potential was improved in ^S] to C'^^^^'^,e > 0, 
and to Dini continuous in ^U]. Recently in ^1] a method was given 
for reconstructing the magnetic field and the electric potential under 
some regularity assumptions on the magnetic potential. 

All of the above mentioned results rely on constructing complex 
geometrical optics solutions, with a linear phase, for the magnetic 
Schrodinger equation. We also mention that the inverse boundary value 
problem is closely related to the inverse scattering problem at a fixed 
energy for the magnetic Schrodinger operator. The latter was studied 
under various regularity assumptions on the magnetic and electrical po- 
tentials in P, for small compactly supported magnetic potential and 
compactly supported electric potential. This result was extended in 
[3] for exponentially decaying magnetic and electric potentials with no 
smallness assumption. 

In this paper we extend the main result of [7^ to the case of the 
magnetic Schrodinger equation. We state the precise results below. 



Let Xf) G R'*\ch(fi) (where ch(f2) denotes the convex hull of Q), we 
define the front and back sides of dfl with respect to Xq by 

F{xo) = {x edQ:{x- Xq) ■ u{x) < 0} 
B{xo) = {x e dQ : {x - Xq) ■ u{x) > 0}. 

Theorem 1.1. Let Q be a simply connected open bounded set with C°° 
boundary in R", n > 3, let Ai, A2 be two real vector fields on Q 
and qi, q2 be two bounded potentials on Q such that the assumption^ 
is satisfied. 

^Here A is viewed as the 1-form J2]=i Ajdxj. 
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Let xq E R"\ch(i7), suppose that the Dirichlet to Neumann maps 
related to the operators CAi,qi and CA^m coincide on part oj the bound- 
ary near xq in the sense that there exists a neighborhood F of the front 
side of dQ with respect to Xq such that 

(1.4) ^fA„,J{x)=AfA„,J{x) \/xeF, WfenHon), 

then if Ai and A2 are viewed as 1-forms 

dAi = dA2 and qi = q2- 

Remark 1.2. We only use the simple connectedness of the set Q to 
deduce that the two magnetic potentials differ from a gradient, from 
the equality dAi = dA2. If we already know that Ai — A2 = V\l/, we 
don't need the fact that Q is simply connected, in particular, theorem 
II. II contains theorem 1.2 of 0. 

Nevertheless theorem 1.1 in [7] improves on this result by restricting 
the Dirichlet-to-Neuman maps to a space of functions on the boundary 
with support in a small neighborhood of the back side B{xo). We have 
left the corresponding result in the magnetic case open. 

As in we make the following definition of a strongly star shaped 
domain. 

Definition 1.3. An open set Q with smooth boundary is said to be 
strongly star shaped with respect to xi G dQ if every line through xi 
which is not contained in the tangent hyperplane cuts the boundary at 
precisely two distinct points Xi and X2 with transversal intersection at 

X2. 

With this definition, theorem 11.11 implies the following corollary 

Corollary 1.4. Under the assumptions on Q, the magnetic potentials 
Ai, A2 and the electric potentials qi, q2 of theorem M . li let xi G dVt be 
a point of the boundary such that the tangent hyperplane of dQ at xi 
only intersects dQ at Xi and such that Q is strongly star shaped with 
respect to Xi . 

Suppose that there exists a neighborhood F of xi in dQ such that 
di.^D holds then 

dAi = dA2 and qi = q2- 

We proceed as in [7^ by constructing some complex geometrical op- 
tics solutions using a Carleman estimate. The construction of these 
solutions is fairly similar to those presented in the latter paper, except 
for the changes due to the presence of the magnetic potential. However, 
the part concerned with the recovery of the potential and the magnetic 
field is new. 
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The plan of this article is as follows. 
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2. Carleman estimate 

Our first step is to construct solutions of the magnetic Schrodinger 
equation CA,qU = of the form 

(2.1) h) = e^^'^+'^\a{x) + hr{x, h)) 

(where ip and ip are real functions) by use of the complex geometrical 
optics method: of course, ip will be sought as solutions of re- 

spectively an eikonal equation and a transport equation. In order to 
be able to go from an approximate solution to an exact solution, one 
wants the conjugated operator 

to be locally solvable in a semi-classical sense, which means its principal 
symbol^ 

(2.2) p^(a;,^) = ^2_^y^)2^2^V<^.^ 

to satisfy Hormander's condition 

{Rej9<^, Imp^} < when = 0. 

Since we furthermore want to obtain solutions (|2.ip for both the 
phases cp and —ip, we will consider phases satisfying the condition 

(2.3) {Rep^, Imp^} = when = 0. 

Definition 2.1. A real smooth function ip on an open set Vt is said 
to he a limiting Carleman weight if it has non-vanishing gradient on 

^Here and throughout this article, we are using the semi-classical convention. 
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and if the symbol satisfies the condition ^2. on T*{VL). This is 

equivalent to say that 

((231) {ip"Vip, Vip) + {ip"i, = when = (Vy?)^ and Vy? ■ ^ = 0. 

The appropriate tool in deducing local solvability for the conjugated 
operator and in proving that the geometrical optics method is effective 
(meaning that indeed one gains one power of h in the former asymp- 
totics) is a Carleman estimate. The goal of this section is to prove such 
an estimate. 

In this section, f2 is as in the introduction and VL will denote an open 
set 3) We will use the following notations 

{u\v) = I u{x)v{x) dx, {u\v)qq = / u{x)v{x) da{x), 
Jn Jan 

and = -y/ {u\u) denotes the norm on fi. We say that the estimate 

Fiu,h)<Giu,h) 

holds for all M G X (where X is a function space, such as L'^{Q)) and for 
h small if there exist constants C > and ho > (possibly depending 
on q and A) such that for all < h < ho and for all m G X, the 
inequality F{u, h) < CG{u, h) is satisfied. 

We will make extensive use of the Green formula for the magnetic 
Schrodinger operator Ca^q-, which for sake of convenience, we state as 
a lemma. 

Lemma 2.2. Let A he a real vector field on VL and q G L'^{VL) then 
we have the magnetic Green formula 

(2.4) {CA,qu\v)^ - {u\CA,qv)n = 

{u\{d^ + iu ■ A)v)g^ - {{d^ + iv ■ A)u\v) 
for all M, f G H^iyi) such that Am, Af G L'^iVl). 
Proof. Integrating by parts, we have 

(2.5) {CA,cu\v)n = iVu\Vv)n + iDu\Av)n + {Au\Dv)n 

+ {{q + A^)u\v)^ - {{d^ + %v ■ A)u\v)q^ 

and by permuting u and f , replacing q by g, and taking the complex 
conjugate of the former, we get 

{LA,ciu\v)n = iVu\Vv)n + {Du,Av)n + iAu\Dv)n 

+ ((g + A^)u\v)^ - {u\id, + %v ■ A)v)q^. 

Subtracting the former to ()2.5|) . we end up with ()2.4|) . □ 
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If 99 is a limiting Carleman weight, we define 

dn± = {x Edn-. ±du^ > 0}. 

Proposition 2.3. Let ip be a limiting Carleman weight on Q, let 
A be a vector field on Cl and q G L°°(i7), the Carleman estimate 

(2.6) - h{dy(p e^ duu\e^ dyu)QQ,^ + ||e^n|p + ||e^WM|p 

< h Weh CA,qu\\^ + h{dyLp eh dyu\eh dyu)dQ,+ 
holds for all u G C°°{Vt) fl -f^o(^) ^'^'^ h small. 

In particular, when u G C^(fi), we have the Carleman estimate 
(12.61 ) + /i||e^Vu|| < h\\eh CA,qu\\. 

Proof. Taking v = e^u, it is equivalent to prove the following a priori 
estimate 

(2.7) - h{d,^dMduv)an- + ll^f + WhVvW" 

< ^\\{ehh^ CA,qe'h)vf + h{d^Lpd^v\d^v)Qn^ 

since ||/ie^VM|| < ||f || + ||/iVf || and {e^dyu)\d^ = d^v\gQ. Conjugating 
the magnetic Schrodinger operator by the exponential weight gives rise 
to the following operator 

(2.8) e^h'^CA,qe~T^ = P + iQ + R + h^{q + A^) 
where P and Q are the self-adjoint operators 

P = h^D^ - (Vy?)^, 
Q = Vip ■ HD + HD ■ Vv9, 
and R = h{A- HD + HD ■ A) + 2ihA ■ V<^. 

Our first remark concerns the fact that we may neglect the term h'^{q + 
A"^) since the right-hand side of (|2.7|) may be perturbed by a term 
bounded by |p, which may be absorbed by the left-hand side if h 
is small enough. Omitting the term q + A"^ gives rise to such an error. 
Hence we will prove the a priori estimate for the operator P + iQ + R. 
The same is not true of the term R because errors of order ||f |p + 
ll/iVf IP may not be absorbed into the left hand-side. 

Note that if p and q denote the principal symbol respectively of P 
and Q, the fact that is a limiting Carleman weight means that 

{P) Q'} = when p + iq = 0. 

This condition is not enough to obtain an a priori estimate for P + iQ, 
one needs to have a positive Poisson bracket. Our first step is to remedy 
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to this by using a classical convexity argument. Consider the modified 
Carleman weight 

where e is a suitable small parameter to be chosen independent of h, 
and denote by p and q the corresponding symbols, and by P,Q,R the 
corresponding operators, when ip has been replaced by (p. Then, we 
have 

= (1 + ^(p)V(p, <p" = (1 + ^^)^" + jV^ (g) Wif 

therefore when = (V<^)^ and V<^ ■ ^ = 0, we have 

(2.9) {p, q} = 4(^"e, + 4(^" V^, V^) 

since (p is a limiting Carleman weight. Furthermore, if we restrict 
ourselves to the hyperplane Vx orthogonal to V^^, we get 

4-h h 

{p, ■)\y^ = -(1 + -^nv^r + - (v^)') 

with a{x) = Ah{y(pY /e — 4((^"V<^, V(f)/(y(py, and since this bracket 
is a quadratic polynomial with no linear part, this implies that there 
exists a linear form b{x, ^) in ^ such that 

{P, <?} = y (1 + ^^fiVip)^ + a{x)p + i)q. 

This computation implies on the operator level that 

~ ~ 4/)2 h h ~ ~ 

i[P, Q] = —{1 + -^fiy^Y + -(aP + Pa) 

where the first order differential operator 6"" is the semi-classical Weyl 
quantization^ of h. In fact, the positivity of the bracket ()2.9|) essentially 
induces the positivity of the commutator i\P , Q\ 

(2.10) 2([P, Q]v\v) = —(1 + -^fWiV^fvr 



e e 



>o 

,3/ 



+ hRe{aPv\v) + hRe{Qv\h'" v) + h\cv\v) 



^The absence of /i^ term is due to the use of the Weyl quantization. 
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(recall that v\qq = which explains why there are no boundary terms). 
The former fact will be enough to obtain the a priori estimate on 
P + iQ. 

Our last observation is that 

llWwf = {Pv\v) + \\^/{Vip)^vf 

leading to 

(2.11) \\hVvf <\\Pvf + \\vf. 
Now, we turn to the proof of the estimate. We have 

\\iP + iQ)vf = \\Pvf + WQvf + i{Qv\Pv) - i{Pv\Qv) 

and the magnetic Green formula (j2.4j) (used in the straightforward 
case with no potential P = /i^>Co _(Vip)2/h2)5 together with the fact that 
v\on = 0, gives 

iQv\Pv) = iPQv\v) - h^iQv\d,v)an 

= {PQv\v) + 2ih^{du(pdyv\dyv)dn 

and similarly, since Q is first order, we get 

{Pv.Qv) = {QPv,v). 

Therefore we have 

\\iP + tQ)vf = \\Pvf + WQvf + ii[P,Q]v\v) 
- 2h^{dyv\du(pdyv)dn 

and using (j2.1U|) . we get 

\\{P + iQ)v\\'' + 2h\d,v\du^d,v)9n > \\Pvf + \\Qvf+ 

+ ll^f - {Ch\\Pv\\ \\v\\ + Ch\\Qv\\ \\hVv\\) 
£ ^ V ' 

<|ll^''P+|||Qi;||2+^(||t,||2+l|Wt>||2) 
which combined with ()2.11|) . gives when e is small enough 

\\{P + iQ)v\\'' + 2h\dMd.(pd,v)a^ 

> (1 - Oie'^h^Wvf + WQvW + ^(||t;f + ||Wi;f ). 
Thus, taking h and e small enough 

(2.12) \\{P + iQ)v\\'' + 2h\d,v\d,^d,v)an > —{\\vf + \\hVvf). 
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The last part^ of the proof is concerned with the additional term Rv 
due to the magnetic potential; from the former inequality we deduce 

II (P + iQ + R)vf + h\dMd.(pd,v)an 



> —I 



V 



|2 



\hVvf -0{e)\\h-^Rvf) 



and using the fact that \\h ^RvW^ < ||f |p + ||/iVf |p, we obtain 

\\{P + tQ + R)vf + h%d,v\d,ipd,v)9n > —{\\vf + \\hVvf) 

if e is chosen small enough. Finally, with v = e^e^u, we get 

1 

7^ 



I ^ £ ii9 II T ^ i£— i£ ii9 ^ II i£— iE , 9 ^ ii9 

|e2eehn|| + ||We2£e'»n|| < —\\e2ee'^hLA,qu\\ 



+ /i(e^= d^{ehu)\d^ipd^{ehu))g^ 
this gives the Carleman estimate ()2.6|) since 

l<e'"'-<C. 1<|^ = 1 + ^<| 

2 Ou^p e 2 

on for all h small enough. □ 

We denote by HI^^{VL) the semi-classical Sobolev space of order 1 on 
VL with associated norm 

= ll^Vwf + ||nf 

and by H^^-^iTi^) the semi-classical Sobolev space on R" with associated 
norm 



\u 



Changing into —(p, we may rewrite the Carleman estimate in the 
following convenient way 

(2.13) Vh\\^d^e~'^dM\ 

< h\\e'hCj^^^u\\ + Vh\\^-dyipe~hdyu\\L^(^da_). 

By regularization, this estimate is still valid for u G H'^{Q) fl Hq{Q). A 
similar Carleman estimate gives the following solvability result: 



^This is the main difference with respect to the proof of the Carleman estimate 
in 13. 
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Proposition 2.4. Let (f be a limiting Carleman weight on Vt, let A he 
a vector field on VL and q G There exists ho such that for all 

< h < ho and for all w G there exists u G H^{VL) such that 

h^CA,q{e^u) = e^w and h\\u\\HijQ) < \\w\\. 

Proof. We need the following Carleman estimate 

(2.14) ||t;|| < /i||e-^£^,,e^i;||^-i(i,„), Vt; G C^iil). 

scl -' 

Let Q <£ Cl G Q, assume that we have extended A to a vector field 
on Q and g to a L°° function on Cl. Let x ^ ^o^i^) equal 1 on Q. With 
the notations used in the proof of proposition 12.31 we have 

{hD)-\P + iQ){hD) =P + iQ + hRi 

where Ri is a semi-classical pseudo-differential operator of order 1, 
therefore from estimate ()2.12|) we deduce 

||(P + ^g)(/.D)^;||5,-. >^(||^;||^, (j,„)-0(e)||i?it;f) 

scl ^ ^ scl ^ -' 

h"^ 

> II l|2 

~ _ Irll/fi, (R") 

for any v G C^(f2), if h and e are small enough. Besides, we have 

\\{R^h^{q^ A^))v\\^-J^^^^^<h\\v\\ 
therefore if e is small enough, we have 

11(^^)^11 ^ ^l|e~^>CA,<^e^(/^,L')^;||^-l/I^„^ 

scl ^ ) 

for any v G Hence if u G C^(J^), taking v = x{hD)~^u G 

C°°{Cl) in the former estimate, and using the fact that 



\\{i-x){hD)-\\\H:^^ = o{h^)\\u\\ 

we obtain 

ll^ll ^ H^^"^ '^A,qe'i^u\\fj-i(^„), \/u G C~(fi). 

scl ^ 

This gives the Carleman estimate ()2.14|) since e'^^^ = e'^^^^e'^^^. Classi- 
cal arguments involving the Hahn-Banach theorem give the solvability 
result. □ 
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3. Construction of solutions by complex geometrical 

OPTICS 

The goal of this section is to construct solutions of the magnetic 
Schrodinger equation of the form ()2.1|) . To do so we take to he a 
solution of the eikonal equation 

p{x, V'?/'(x)) + iq{x, V^J^x)) = 

such solutions exist since {p, q} = when p = q = 0. More precisely, 
the eikonal equation reads 

(3.1) (V^)2 = (V(^)2^ V<^-VV' = 0. 

In fact, in the remainder of this article, we fix the limiting Carleman 
weight to be 

(3.2) ip{x) = ^\og{x - xof. 

For such a choice of ip, the second part of the eikonal equation is 
merely the fact that ip is a. function of the angular variable {x — xo)/\x — 
Xo\ and we can actually give an explicit solution of the eikonal equation 

(3.3) ip{x) = arctan ^ 



2 a/(x - Xo)2 -{u- {x- Xo)y 

X — Xq 



\X — Xo\ 

where u G S*""^. Let us be more precise about the set where u may 
vary, keeping in mind that we want this function to be smooth — in 
particular, we have to ensure that u; 7^ (x — Xq)/\x — xo\ whenever 
X eCl. 

For that purpose, let tq > be large enough so that Q C B{xo, tq), 
let H denote a hyperplane separating xq and ch(fi), and the open 
half space containing Q (and therefore xq ^ H^), we set 

r = G 5"-^ : Xo + rod G H+} 

and r the image of F under the antipodal apphcation. Let ojo G 
Sn-i\{Y u f ) and Fq be a neighborhood of uq in S"-i\(F U f ), then 
the distance F x Fq 9 {0,uj) —>■ dsn-i{6,uj) is a function. Moreover, 
fl (Z fl = xo + R+F hence we have {x — xo)/\x — xq | G F for all x in the 
open neighborhood Cl of 0, thus iJj depends smoothly on the variables 
(x, uj) on Q X Fq. 

Remark 3.1. Suppose that Xq = and u = (1, 0, . . . , 0), which we can 
always assume by doing a translation and a rotation. Notice that by 
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considering the complex variable z = Xi+i\x'\ & C (with x = {xi,x') G 
R X R""^), we have 

(p = log \ z\ = Re log z, ip = arctan = Im log z 

2 Imz 

when Imz > (note that ip = arctan(Im2;/ Kez) on the first quadrant 
Rez > 0, Imz > 0) hence ip + iip = log z. 

With such (fi and ip, we have 

h^JZA,qei^''^'^^ = er^^^'^'^ [h{D + A) ■ (W - iV^) 

+ hiVtlj - iVip) ■{D + A) + h^CA,q) 

thus we will have 

h^CA,,{e^^^+''^^a) = 0(/i2)e^(^+*'^) 
if a is a solution of the first transport equation, given by 

{{D + A) ■ (VV^ - iVip) + (VV^ - iVip) ■ {D + A))a = 0. 
We write the latter vector field equation 

(3.4) (V^ - iVif) ■ Da + (Vip - iVip) ■ Aa + —(Aip - iAif)a = 0. 

2i 

We seek a under exponential form a = e*, which means finding $ 
solution of 

(3.5) (V(/? + iV^) ■ V$ + + iVi)) ■ A + ]^A{ip + i^)) = 0. 

on Q. The function $ has regularity since the magnetic potential A 
is C\ 

Remark 3.2. Considering the complex variable z = xi + i\x'\ as in 
remark [3.11 with (f + tip = log 2;, we may seek $ as a solution of the 
following Cauchy-Riemann equation in the z variable 

a$ (n-2) 1 , 

where = (0, 6) is the unit vector pointing in the direction of the 
r-axis. Indeed, if we denote by (xi, r,9) G R x R+ x S*""^ a choice of 
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cylindrical coordinates on R" and z = Xi + zr, we have 

dlogz d dlogz d 2d 



dxi dxi Or dr z dz 
92 {n-2) d 1 



/ o o [n — Z) o 1 \ 
and A(v, + = + ^ + —^Q-, + -2^s^-^) + 

n-2 d ^ (n-2)i 

= Trlog^: = • 

r or rz 

Remark 3.3. Note that the set of solutions of (j3.4p is invariant under 
the multiplication by a function g satisfying 

{\/(fi + iVip) ■Vg = 0. 

In the setting of remark this condition reads 

oz 

on Q, i.e. g is a holomorphic function of z = xi + i\x'\. 

Having chosen the phase if + and the amplitude e*, we obtain an 
approximate solution of the magnetic Schrodinger equation 



9\ 



(recall that $ is C^) which we can transform into an exact solution 
thanks to proposition 12.41 there exists r(x, h) G H^{Vt) such that 

and MH^^^in) < ||>C^,ge*||. 

We sum up the result of this section in the following lemma. 



Lemma 3.4. Let xq E R"'\chr2, there exists ho > and r such that 
Mn^jn) = and 

u{x, h) = e^(^+^'^)(e*(^) + hr{x, h)) 

is a solution of the equation CA,qU = 0, when h < Hq, and (p is the 
limiting Carleman weight (j.V.^ . ip is given by fl.V. 'J^ and ^ is a solution 
of the Cauchy-Riemann equation (j^'.^p . 

Note that with ip as in ()3.2j) the parts of the boundary dQ± delim- 
ited by the sign of di^(p correspond to the front and back sides of the 
boundary 
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4. Towards recovering the magnetic field 

Let Xq G R"\fi, suppose that the assumptions of theorem 11.11 are 
fulfilled and consider 

= {x G dVl : (x — Xq) ■ i^{x) < 6\x — Xop} D F{xo) 

with e > small enough so that C F, therefore satisfying 

(4.1) J\fA„,J{x)=AfA,,gJ{x), VxGF,, WfeH'HdQ). 

We may assume without loss of generality that the normal compo- 
nents of Ai and A2 are equal on the boundary 

(4.2) Ai-u = A2-iy ondn 

since we can do a gauge transformation in the magnetic potential 

(see (|1.3|) in the introduction) with \E' G 0^(0,) such that ^^I^q = and 
di,'^ is a prescribed function on the boundary^. We extend Ai and 
A2 as compactly supported^ functions in R". 
We consider two geometrical optics solutions 

Uj{x, h) = e^(^^+^^^)(e*^ + hr^ix, h)) 

of the equations Ca^^qiUi = and CA2,q2'^2 = constructed in the 
former section with phases 

(4.3) ip2ix) = -(^i{x) = ip{x) = log \x-Xo\, 
i^2{x) = ipi{x) = ipix) = dsn-1 (-1^ — -H"'^)' 

\\X — Xq\ / 

defined on a neighborhood of (and uj varies in To), and where $1 
and $2 are solutions of the equations 

(4.4) (Vy? - iVij) ■ V$i + i(V(/? - iVij) ■ Ai + ^A{ip - i^p) = 

(V<^ + iV'^) ■ V$2 + i(V<^ + iVtp) ■ A2 + ^A{!f + iip) = 0. 

The remainders vj are bounded independently of h in H^^^. Note that 
it implies the following estimate on Uj 

(4.5) \\e-^u,\\H^^^ = Oil). 



^After the use of a partition of unity and a transfer to {xi > 0}, this is theorem 
1.3.3. in 0. 

^Note that Ai and A2 do not necessarily agree on dfl. 
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By w we denote the solution to the equation 

Ca^qiW = 0, w\dn = U2\dn 

so that AfAi,qi{u2\dn) = {duw)\dQ + iAi ■ i'U2\dn- The assumption dHH) 
means that 

d,y{w — U2){x) =0, Vx G 

(here we use the fact (|4.2j) that the normal components of the magnetic 
potentials coincide on the boundary). Besides, we have 

(4.6) £Ai,gi(w - U2) = -CAi,qiU2 = {CA2,q2 - ^Ai,qi)u2 

= {A2 - Ai) ■ Du2 + D-{A2- Ai)u2 
+ {Al -Aj + q2- qi)u2 

hence we deduce 

(4.7) {CA,,qAw - U2)\Ul) 

= {{A2 - A) ■ Du2\u^) + {u2\iA2 - A,) ■ Dm) 

+ - ((^2 - Al) ■VU2\Ui) +[{Al -A\ + q2- qi)u2\ui). 

^ ' V ' 

=0 

The magnetic Green's formula gives 

(4.8) {CAuqA^ - U2)\Ui) 

= {w - U2\CA^,qiUi) -Udy + iAi ■u){w - U2)\ui)dn 

' V ' 

=0 

= -{d^iw - U2)\ui)dn\F, 
and combining (j4.7p and (j4.8|) . we finally obtain 

(4.9) / du{u2- w)uida{x) = j {Al- A\ + q2 - qi)u2Uidx 

+ / (A2 - Al) ■ {Du2 ui + U2 Vui) dx. 
Jn 

With our choice of = = log |x — Xo|, we have 

F, D F{xo) = on. thus dn\F, C dn+ 

and moreover d^^p > e on dQ\Fe therefore the modulus of the left-hand 
side in ()4.9p is bounded by 

-^\\>/d^e'hd^{u2 - w)\\dn^ x ||e*^ + hriWgn^ 



<l|e*i||af2 + ||r-i 
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which, in virtue of the Carleman estimate (j2.13p . is bounded by a 
constant times 




=0 because of 



In view of ()4.f)|l and of ()4.5j) the former expression is 0{h 2 ). Therefore 
we can conclude that the right-hand side of ()4.9|) is 0{h~2). This 
constitutes an important difference with 0, where the corresponding 
term was 0{h). 

More directly, the first and second right-hand side terms of ()4.9|1 are 
respectively and 0{h~^) as may be seen from (|4.5p . It turns out 
that the information obtained when disregarding the bounded term is 
enough to recover the magnetic field. We multiply ()4.9|) by h and let 
h tend to 0: 

lim / ((^2 - Ai) ■ hDu2 ui + U2 {A2 - Ai) ■ hDui) dx = 0. 
Using the explicit form of the solutions ui and U2, this further means 

(4.10) / {A2 - Ai) ■ (V<^ + zV^/')e*^+*2 dx = 0. 

Adding the complex conjugate of the first line of ()4.4|) to the second 
line, we see that 

(Vy? + iVip) ■ (Vi-i + V$2 + i{A2 - Ai)) + A{ip + iilj) = 

this implies that 

(4.11) {D + A2- Ai) ■ (Vy? + iVV')(e*^+*^) = 0. 

As observed in remark in the expression for U2, we may replace 
e*^ by e'^^g if (7 is a solution of (Vy? + iVip) ■ Vg = 0. Then ()4.10|1 can 
be replaced by 

[ {A2 - Ai) ■ (V<^ + iV^)e*^+*^^ dx = 0. 
Jn 

From equation (|4.1ip . we see that we can replace A2 — Ai by iV in the 
former equality 

(4.12) [ g{x) V ■ (e*^+*^(V<^ + iV^j)) dx = 
Jn 

for all functions g such that (Vf + iVip) ■ V(7 = on ri. 
'''See (5.16) and the subsequent lines in [7]. 
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Remark 4.1. An integration by parts gives 

[ {d^^ + id^i))e^'+^^g da - ! 6*1+*^ ( V(/? + iS/i)) -Vg dx = 
Jan Jn v ■' 

=0 

hence we have / {d,y(p + id^ip)e'^^~^^'^ g da = 0. 
Jan 

5. Moving to the complex plane 

In this section, we follow remark I^TD and choose to work in the cylin- 
drical coordinates. Let us be more precise: suppose that xq = ^ 
ch(f2) and that we have picked u G S'"~"'^\(r U F) with the notations of 
sectional After a rotation, we assume that u = (1, 0, . . . , 0), therefore 
we have 

^] (E C {x G R" : 0}. 
We choose the following cylindrical coordinates 



t = xu r=|x'|>0, e = e(x) = ^es 

\x'\ 



n-2 



By Sard's theorem, the set of critical values of G : t— > 5"""^ is of 
measure 0, therefore the set fle^ = Q~^{9o) = {x E Q : x' = r9o, r > 0} 
is an open set with smooth boundary for almost every in O(f^). The 
result obtained in the former section reads 

/ // 9{x) V, ■ (e*^+*HV,¥P + ^V,^)) r^-^drdtde = 0. 
Js'"-^ J Jug 

and taking g = gi{t, r) ® g2{0) and varying g2 leads to 

(5.1) // (7(t,r)V,- (e*i+*^(V,</^ + 2V,V^))r"-2M = 



for any function g such that (Vxf + i^xi') ■ ^xg = on Qg, and this 
for almost every 6. 

Now we consider the complex variable z = t + ir E C+ = {w G C : 
Imw > 0} and write our results in this setting. Let us recall the results 
of the computations made in remarks 13.11 and 13.21 

2d in — 2)i 

(5.2) V^((^ + #) ■ V, = -— and A,(^ + #) = ^ ^ 

z oz rz 

2d in — 2)i 

and thus ■ o (v^ + i'lp) = -ttt H • 

z oz rz 

Similarly, the functions $j satisfy 

(5.3) |,(<li + $2) = + ^(^1 - A,) . (ei + ler) 
oz [z — z) 2 
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Finally, (j5.1|) reads 

for any g G 1-L{VLq). Replacing the holomorphic function g/ z onVLghy 
g, we can drop the factor 1/z. 

If (7 is a holomorphic function, we have^ 

d{e^'+^''g{z){z~z)^'^dz) 

= ^ ((2 - zr'^e^'+''')g{z)dz A dz 

therefore the Stokes' formula implies 



/ g{z)e^'+'''{z-z)''-^dz. 
Jane 



Together with ()5.4p this gives 



(5.5) / (7(2)e*i+*^ {z - zy-^ dz = Q 

JdUg 

for any g G TC^Qe). 

Lemma 5.1. There exists a non-vanishing holomorphic function F 
on Qg, continuous on Cle, whose restriction to dVLg is equal to {z — 

Proof. We denote f{z) = {z — ^)"'~2e*^"'"*^ and consider the Cauchy 
integral operator 

c{f)i^) = 7^f P^dc, yzec\dng. 

The function C{f) is holomorphic inside and outside VLq and the Plemelj- 
Sokhotski-Privalov formula reads on the boundary 

(5.6) jim CU){z) - ]im C(/)(z) = f{z,), V^o G d^e. 



^The result of this computation is the transcription of the fact that the formal 
adjoint of (Viy9 + iV^A) • V is V • (V(p + iVijj) in the complex setting, where the 
measure is r"~'^dz A dz. 
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The function ( — z)~^ is holomorphic on Qg when z ^ Qg hence 
(|5.5|) imphes that C{f){z) = when z ^ fig. The second hmit in (j5.6j) 
is 0, thus F = C{f) is a holomorphic function on Qg whose restriction 
to the boundary agrees with /. 

It remains to prove that F does not vanish on Qg. This is clear by 
the argument principle since 

vararg F = vararg / = 

dflg d^lg 

and F is holomorphic. □ 
In particular, with the former function, we have on the boundary 

^ [z - z)"-2e*i+*2 = 1, vz G ang 



F{z) 
which implies 

l-i + $2 + \og{z - = log F{z), \fz e OQe 
with log F a holomorphic function on Qg , and therefore 

/ (/(;2)(<li + $2 + log(z-z)"-2)dz = 0. 
Jang ^ ^ 



An application of Stokes' formula gives 

n-2 



z — z 

hence using the equation this implies 



dz A dz = 



(5.7) // g{z){Ai - A2) ■ {ei + ier) dz A dz = 0. 

J JQg 

With g = 1, the former equality reads 

{t + ir){V^(p + iV^ip) ■ {Ai - A2) dtdr = 0. 



Denote by Pg = span(cij, e,.) the plane along the axis directed by 
LO = (1,0,..., 0), and by P^^ the half plane where x ■ > 0, then 
ng = nn{x= (xi, x') e R" : x' = r^, r > 0} = O n Pfj-. Let TTg be the 
projection on Pg and dXg the measure on the plane, then ()5.7|) (with 
g = I) imphes 

/ ng{Ai- A2)dXg = 0, 

J Pgnn 
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for almost every 9 G S*""^, hence for all 9 G S*"^^ by continuity. The 
former may be rephrased under the form 

f e- (lf,(Ai-A2))rfAp = 0, V^gP 

Jxo+P 

for all linear planes P containing uj = (1, 0, . . . , 0). We can also let Xq 
vary in a small neighborhood of ^ ch(i7) and u vary in a neighborhood 
To of (1, 0, ... , 0) on the sphere 5""^ 

Lemma 5.2. Let A be a vector field on Q. If 
(5.8) / ^■AdXp = 0, VeGT,.(P) 

for all planes P such that d{{0, Ci), T(P)) < 6 then dA = on Q. 

The proof of this lemma is based on the following microfocal version 
of Helgason's support theorem. 

Theorem 5.3. Let f G C°(R"), suppose that the Radon transform of 
f satisfies 

nf{H)= [ fd\H = o 

Jh 

for all hyperplanes H in some neighbourhood of a hyperplane Hq then 

N*{Ho) nWFM) = 
where N*{Hq) cenotes the conormal bundle of Hq. 

The proof of this result may be found in (see proposition 1) or in 
[H] (see section 6). We will also need the microfocal version of Holm- 
gren's theorem (see section 1 or JB] section 8.5). 

Theorem 5.4. Let f G £'(R") then we have 

iV(supp/) CWF,(/) 

where A^(supp /) is the normal set of the support of f . 

These two results may be combined to provide a proof of Helgason's 
support theorem (see [2] and jS]). We also refer to the book j|4j for a 
review on Radon transforms. 

Proof of lemma (3^ Let x ^ ^o^d^l < |) ^"^^ Xe = ^ "x(V^) be a 
standard regularization, one has 

[ ^-iXe* InA) dXp = e-" / x(-) ( / ^ ■ AdX^y+p) dy = 

Jp J ^ ^J{-y+p)nn ^ 
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when (i((0, ei), T(P)) < 6 — e. Therefore it suffices to prove the result 
when Q = W and Ae (R"; R") since d{xe * InA) tends to dA as a 
distribution when e tends to 0. 
Our ffist step is to prove that 

(5.9) L*HdA = 

for any subspace H C R" of dimension 3 such that d{{0, ei),T{H)) < S. 
Here lh denotes the injection of H in R". For any plane P G H such 
that (i((0, ei), T(P)) < 6, we have 

[ {dA,^ AT]) dXp = ^( [ r]-AdXp- [ ^-AdXp^ 

when C,,"!] & Tx{H). The space H is of dimension 3 so we can assume 
that either 77 or ^ belongs to Tx{P), thus the former expression is zero 
because of (jS.Sj) and of the fact that whenever 1] G T^^P) 



t=o 



^ ■ AdXp is constant. 

Therefore, if TZh denotes the Radon transform in H, we obtain 

nH{{L*HdA,^A7])){P)=0 

for any plane P C H such that d{{0, Ci), T{P)) < 6 and for any ^, 77 G 
Tx{H). Combining theorems 15.31 and 15.41 we obtain 

N*{P) n A^(supp(4rfA, ^ A r/)) = 

for any plane P C H such that (i((0, Ci), T(P)) < 6. This gives (j5.9j) 
since such a family of planes sweep other Qn H and the support of A 
is on one side of at least one such plane. 
The result ()5.9p implies in particular that 

{dA{x),^Ar]) = 0, Vx G W,y{^,r]) e S"" x |^-ei| < 5 

and therefore dA = by linearity. □ 

6. Recovering the potential 

End of the proof of theorem \l.l\ Applying this lemma, we finally ob- 
tain 

dAi = dA2 on Q 

therefore the difference of the two potentials is a gradient A1—A2 = V^E' 
(recall that Q is simply connected). The identity (j5.7p now reads 

g{z)d2'^{z,9) dz Adz = 
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for any hofomorphic function g on fig and by Stokes' theorem we get 

/ g{z)'${z,e)dz = 0. 

Reasoning as in the beginning of lemma I^TTl there exists a hofomorphic 
function \& G TC{Qg) such that ^\dng = "^Idfie- Now \l/ is real- valued, 
and since \E' is real-valued on dQg and harmonic, it is real-valued ev- 
erywhere. The only real-valued hofomorphic functions are the constant 
ones, so \l/ and hence is constant on dQe. Varying 6 and also slightly 
Xo and uj, we get that \E' is constant on the boundary dQ. We may 
assume then that "^lan = 0. 

By a gauge transformation, we may assume that ^ = 0, thus Ai = 
Ai- We could almost directly apply the result in [7j to recover the 
identity of the two potentials qi = q2, if it were not for the presence 
of the two magnetic potentials in the equations. Instead we go back 
to the limit induced by (|4.9p . The second right-hand side term is now 
zero. The left-hand side is now 0{\/h) since the 0{h^^) term in ()4.fij) is 
zero and we can reproduce the arguments given after ()4.9|) . Therefore 
we obtain 

(6.1) lim / (g2 — qi)u2Ui dx = 



thus 



/ {q2 - gi)e*i+*^ dx = 0. 
Jn 



As observed in remark 1X51 we may replace e*^ by e'^'^g if is a solution 
of (Vf + iVip) ■ V(7 = 0. Then the former can be replaced by 



(^2 — Q'l) 6*^^*^(7(0;) dx = 0. 



In 

Moving to the complex plane, as in section |3[ this reads 



{q2 - qi)g{z)e'^'+^^{z - z)""-^ dz A dz = 

ng 



for any hofomorphic function g onQg. But the transport equation (|5.3|) 
now reads 

|-((z-z)"-2e*^+*^) =0 
oz 

therefore if we take g = {z — ^)^"+^e~*^~*2, we obtain 



// q{t,r,e)dtdr = 
J J no 
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with q = qi — q2- As in section El varying Xq and u, this may be 
interpreted as 



Jp 

for any plane P such that (i((0, ei), T(P)) < 6. This imphes that for 
any subspace H of dimension 3 such that d{{0, ei),T{H)) < S we have 



for any plane such that d{{0, ei), T(P)) < 6. Regularizing and applying 
theorems 15.31 and 15.41 as in the proof of lemma 15.21 we get 



and therefore l^g = on if leading to g = on fi. This ends the proof 
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